Quasi-Eigenstate Evolution in Open Chaotic Billiards 
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We experimentally studied evolution of quasi-eigenmodes as classical dynamics undergoing a transi- 
tion from being regular to chaotic in open quantum billiards. In a deformation-variable microcavity 
we traced all high-Q cavity modes in a wide range of frequency as the cavity deformation increased. 
By employing an internal parameter we were able to obtain a mode-dynamics diagram at a given 
deformation, showing avoided crossings between different mode groups, and could directly observe 
the coupling strengths induced by ray chaos among encountering modes. We also show that the 
observed mode-dynamics diagrams reflect the underlying classical ray dynamics in the phase space. 
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Quantum manifestation in a classically chaotic system 
has become an important issue in atomic, nano, meso- 
scopic physics, etc., due to its fundamental importance in 
quantum mechanics and applications to practical quan- 
tum/wave systems Most of early works have focused 
on statistical analysis of eigenvalues and eigcnfunctions 
and comparison with the random matrix theory, e.g., the 
transition from Poisson to Wigner distribution of level 
spacings during a transition to chaos, providing an av- 
eraged view on mode dynamics Experimental ver- 
ifications of the statistics have been performed mainly 
in closed microwave cavities Dynamical tunneling or 
coupling between regular and chaotic modes has recently 
been observed for a mixed phase space specially tailored 
for this purpose 0]. 

In open quantum systems, each quasi-eigenmode has a 
linewidth, and thereby changes the mode dynamics sig- 
nificantly. Trapped modes were observed showing high Q 
even with increasing coupling strength to open channels 
in microwave cavities Q , and crossing and avoided cross- 
ing (AC) of cavity modes were reported near an excep- 
tional point formed by two coupled microwave cavities 
0. We note, however, that the previous experimental 
works in microwave cavities and other systems neither 
realized an optimal system showing a continuous chaotic 
transition from being regular to chaotic nor provide ob- 
servations direct enough to tell the variation of statistics. 

In this paper, we have experimentally observed, for 
the first time, the evolution of quasi-eigenmode dynam- 
ics in a generic open nonintegrable system when classical 
dynamics undergoes a transition from being regular to 
fully chaotic. In a dielectric deformation- variable chaotic 
optical microcavity (COM) we traced all high-Q cavity 
modes in a wide range of frequency as the cavity deforma- 
tion increases. By introducing an additional parameter 
orthogonal to the cavity deformation, we could explicitly 
observe mode-mode dynamics under the chaotic transi- 



tion and measure various mode-mode coupling constants 
which can be associated with the underlying classical ray 
dynamics in phase space. We believe our data would 
be a valuable asset for future formulation of a currcntly- 
noncxisting semiclassical theory for coupling strengths 
between modes in a mixed phase space. 

Our experiment was performed in a two-dimensional 
COM made of a liquid jet column]^ of ethanol (refractive 
index m= 1.361) doped with either Rhodamine B dye at a 
concentration of lO^^mol/cm'^ or Rhodamine 6G dye at 
a concentration of lO^^mol/cm'^, depending on the wave- 
length region of interest. Its boundary is approximated 
by r(0) ~ a(l -|- rjCOs24> + erf cos4(/)) in the polar coor- 
dinates with a ~14.9±0.lAim and e = 0.42 ± 0.05 0,0]. 
The deformation parameter 77 can be continuously var- 
ied from 0% to 26%. The size parameter defined as mka 
with k ~ 2ti/\ and wavelength A '^600 nm is about 200, 
thus comprising the short wavelength limit. We mea- 
sured cavity-modified fluorescence (CMF) and/or lasing 
spectra by using the method described in Refs. [ol. [lo|. 

Let us first examine a part of spectrum obtained for 
77=18.7% as shown in Fig. [T]Ja), where each peak cor- 
responds to a cavity mode or a quasi-eigenmode of the 
deformed cavity. The spectrum in Fig. [TJa) consists of 
five different mode sequences. Modes in each sequence, 
marked by vertical ticks below the spectrum, are sep- 
arated by a well-defined interval Aj/ similar to regular 
modes in a symmetric cavity. This is because all of these 
modes are far apart accidently in this frequency region 
and thus any possible interactions among them can be 
neglected. We call them uncoupled. In this limited range 
of frequency we can then label these uncoupled mode se- 
quences by mode order lo{=l, 2, . . . , 5) in the increasing 
order of their FSRs {Ai/i < A1/2 < ■ • ■ < A1/5) in analogy 
to the radial quantum number for a circular cavity [l0|. 

Outside the frequency range of Fig. [TJa), however, 
some of the modes from different mode sequences would 
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FIG. 1: (Color on line) (a) For a COM with 77=18.7%, five 
uncoupled mode groups are identified, (b) Several AGs are 
observed as we vary rj when modes are followed adiabatically. 
Inset: the spectrum measured with a spectrometer with a 
higher resolution (~0.012 THz), resolving AC of two adjacent 
modes. 

get very close because of their different Av's and they 
would interact and repel each other due to the coupling 
introduced by ray chaos as to be seen later. Even in this 
case, we can extend our labeling over the entire spectral 
range of measurement (420 THz - 530 THz) by employ- 
ing the conventional assumption of adiabatic change in 
Av of a given mode sequence from one FSR to another. 
In order to distinguish this modc-scqucncc label from the 
mode order defined above, we use a different notation /, 
called mode label, such that / coincide with Iq only in the 
above limited region of frequency. 

CMF and lasing spectra similar to that in Fig. [TJa) 
have been measured for various rj values from 10% to 
23% and all of the observed modes are labeled by the 
convention explained above. A part of the results are 
shown in Fig. [TJb) , where we can see some of encounter- 
ing modes {1=2 mode and 1=4 mode, 1=1 mode and 1=2 
mode) undergo AGs as the cavity deformation is varied. 

In order to investigate mode dynamical properties at 
a fixed deformation, we now introduce an internal pa- 
rameter n indexing the recurring modes in a given mode 
sequence. The usefulness of n is obvious when we con- 
sider the quasi-eigenmodes of a deformed cavity, obtained 



by diagonalizing a two-state effective Hamiltonian matrix 
(with h=l), 

^ ' [ CM Mv)-^iM\' ^ ' 

where Vpi^^^ and 7p(q) are frequencies and decay rates 
of two uncoupled states with different mode orders, re- 
spectively, and Cqp is the internal coupling induced by 
cavity deformation. This coupling is taken to be real 
because it arises mainly from internal ray dynamics in 
our experiment as to be shown later. AC along rj as 
shown in Fig. 1(b) then takes place at r^o satisfying 
^p{n,m) = ^q{n,m) for a given n if Cpq > |7p - 7g|/2, 
a criterion for AC. Now we consider the variation of n 
at a given rj'{^ rjo) instead. Since states with different 
mode orders have different Ai^'s, there exists some tiq 
satisfying Vp{no,r]') ~ Vg{no,ri'), for which an AC can 
take place. We assume that both decay rate 7p(g) and 
coupling strength Cpq are independent of n since their 
dependence on frequency is not substantial in the fre- 
quency range studied. 

Mode-dynamics diagrams in Fig. [2] are based on this 
idea of scanning n. We first define reference frequencies 
as the resonance frequencies of Zo=3 whispering-gallery 
modes in a circular cavity whose round trip length is the 
same as that of the COM under investigation. These ref- 
erence frequencies are shown as equally-spaced vertical 
ticks marked as 'ref in Fig. [IJa). We then measure the 
relative frequencies of the observed quasi-eigenmodes cor- 
responding to n with respect to the reference frequency 
of the same n for a given r], and plot these relative fre- 
quencies as a function of the reference frequency corre- 
sponding to n. Mode-dynamics can be analyzed more ef- 
fectively in a mode-dynamics diagram than in Fig. [TJb) 
since we can then associate the observed mode dynam- 
ics to the relevant phase-space structure for intracavity 
ray dynamics, the so-called Poincarc surface of section 
(PSOS), for a given ry. 

Note in Figs. [2I^b)-(d) that when these quasi- 
eigenmodes are far apart they follow straight lines called 
diabatic transition lines [llj even in the presence of the 
internal coupling C [the case of Fig. 1(a)]. By shifting 
the internal parameter n, we can bring any two quasi- 
eigenmodes get close and make the internal coupling 
come into play. In this case, the quasi-eigenvalues devi- 
ate from the diabatic lines significantly, exhibiting ACs. 
Note also that the mode order Iq is associated with the 
uncoupled states located on the diabatic lines (straight 
lines in Fig. 2), while the mode index / is associated with 
quasi-eigenmodes on adiabatic lines (exhibiting ACs in 
Fig. 2). The shorthand notation such as 1^ in Fig. 
2 is based on this idea. Furthermore, by comparing Fig. 
^a) in the case of circle with Figs.[2ljb)-l2{d) for deformed 
cavities, we can recognize that the modes on the ^oth di- 
abatic line must have evolved from the WGM's of radial 
quantum number Iq of a circular cavity. 
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FIG. 2: (Color online) (a) Calculated relative frequencies of 
quasi-eigenmodes with radial mode order Zo=l, 2, ■ ■ ■, 8 for 
a circular cavity (77=0%) with the mode frequency of Zo=3 
as a reference, (b) Observed relative frequencies of quasi- 
eigenmodes for ?7=14.3%. Mode frequencies more or less fol- 
low diabatic lines (dotted straight lines) except for ACs with 
very small splittings. We employ shorthand notation Z'" as 
explained in the text, (c) The same for »7=18.7%. More pro- 
nounced ACs with decay-rate exchange as well as ordinary 
crossings are observed. The diameter of the circle drawn on 
each data point represents the half linewidth of the corre- 
sponding mode in THz. Red circle indicates spectrometer 
resolution, 70 ~ 0.05 THz. (d) The case of 77=22.3%. The 
splittings are much more larger than those of (c) and the mode 
frequencies deviate greatly from the diabatic lines. 



The diameter of the circle drawn on each data point 
in Fig.[2{c) represents the half linewidth in THz, directly 
observed with a spectrometer. It is reassuring to see 
that the linewidth well before and well after an avoided 
crossing is continuous along the diabatic transition line, 
which is a general property of avoided crossing On 
the other hand, in the region where avoided crossings 
occur, the linewidth is an intermediate value of those 
well before and well after the avoided crossing. 

Another important factor to consider in Fig. [5] is the 
parity of mode. Only modes with the same parity can in- 
teract with each other. In the frequency range of ly ~500 
THz and < i// — i/rcf <2.25 THz, uncoupled states of 
Zo=3 and 5 have a parity different from that of /o=li 
2, and 4 states of the same n. This feature has been 
confirmed by mode calculations by boundary element 
method [l^, This is why quasi-eigcnmodcs origi- 

nating from lo=l, 2 and 4 states avoid each other there 



FIG. 3: (Color online) (a) Decay rates 7p of p=l,2,4 uncou- 
pled states, (b) Coupling strength Cpq between p, q(=l,2,4) 
uncoupled states, restored from the observed sizes of AC and 
the decay rates of modes, assuming three-mode coupling. Sys- 
tematic error in determining exact deformation is indicated 
with a horizontal error bar in both (a) and (b). (c) PSOS' for 
rj— 17% and (d) r]=19%. Large (purple) circular and (orange) 
square dots represent the classical trajectories that Iq — I and 
2 modes would correspond to, respectively, whereas that of 
lo=4 mode is embed in the chaotic sea below. Inset: Birkhoff 
coordinates used in PSOS'. 



and why Zo=l and 2 states cross the lo=3 state near 
(j/, i^i - t'ref)'-(510, 0)(THz) in Figs.[2i;b)-ii;d). However, 
the same lo=l state and another lo~3 state displaced by 
one FSR result in quasi-eigenstates undergoing an AC 
near (430, 2.3) (THz) since any state with its mode num- 
ber shifted by one (n — > n ± 1) would have its parity 
changed to the other parity [l^ . 

By the same reason one may expect that the same Iq=1 
state and another la=5 state with an one-less mode num- 
ber would result in an AC near (525, -0.25) (THz) in Figs. 
He), (500, 0.25)(THz) in Figs.ElJb), and (540, 1.4)(THz) 
in Figs. El^d). However, they all exhibit a crossing in- 
stead. It is because C15 < I75 — 7i|/2, not satisfying the 
criterion for AC. This example demonstrates that open- 
ness can suppress the AC in the present internal coupling 
case. From this openness effect, we can expect that the 
level spacing distribution would show a delayed transition 
from Poisson to Wigner-like distribution in the chaotic 
transition. 

From the observed gaps of AC and the associated de- 
cay rates of corresponding uncoupled states [Fig. [3]Ja)], 
we can finally reconstruct the internal coupling strength 
C{r]) between encountering modes as shown in Fig. 
^h). In the present case, all three ^o=l,2,4 modes of 
the same parity are coupled to each other since their 
mode frequencies are not much separated in the region 
of interaction. The reconstructed coupling strengths. 
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FIG. 4: (Color online) (a) Calculated relative frequencies of 
1=1, 2, 3, 4, 5 modes with respect to the reference frequency, 
(b) Hisimi plots of the modes marked by arrows near 460 THz 
where 1=2, 4 modes undergo an AC. 



Ci2, C24, (^14, summarized in Fig. [3]Jb) are obtained by 
diagonalizing a three-mode non-Hcrmitian symmetric 
Hamiltonian. a straightforward extension of Eq. (1). It 
is the first time to directly measure mode-mode coupling 
constants in an open chaotic billiard of generic noninte- 
grable shape. In Fig. [3Jb) these couplings are shown to 
increase as the degree of deformation increases. 

Unfortunately, there is no known semiclassical theory 
for enabling us to calculate the observed coupling con- 
stants. At best, they can be understood qualitatively in 
terms of classical ray dynamics in phase space. Follow- 
ing this standard practice we plot PSOS in Figs. [Sj^c) 
and[3jd), for r/=17% and 19%, respectively, by using the 
Birkhoff coordinates with the polar angle and x the 
incident angle in ray tracing analysis 0] . Large (purple) 
circular and (orange) square dots represent the classi- 
cal trajectories that /o=l and 2 modes would correspond 
to, respectively, whereas that of Zo=4 mode is embedded 
in the chaotic sea for the shown degrees of deformation. 
These trajectories are inferred from phase-space distri- 
butions or Husimi plots of the these modes [lO|. When 
T] > 18%, as shown in Fig. [3Kd), the classical trajectory 
associated with lo=2 mode no longer lie on the main inte- 
grable region, separated from the chaotic sea by unbro- 
ken Kolmogorov-Arnold-Moser (KAM) curve, as it did 
in Fig. [3jc) for r/=17%, but lie on islands surrounded by 
chaotic sea, and thus chaotic diffusion starts to play an 
important role for the increased coupling C24 between 
lo=2 and 4 modes as shown in Fig. [Sja). The broken 
KAM curve is also responsible for the increased coupling 
Ci4 between la=l and 4 modes. 



The observed mode-dynamics diagrams have also been 
reproduced by numerical calculations based on the 
boundary element method 12, 3] applied for the same 
shape and size of the cavity as in the experiment. The 
eigenvalues and associated Husimi distributions calcu- 
lated for 77=0.19 are shown in Fig. [H where we confirm 
that the encountering quasi-eigenmodes exchange their 
mode distributions upon avoided crossing [(1)^(6) and 
(4)^(3)] and at the closest encounter the resulting modes 
[(2) and (5)] are linear superpositions of the modes well 
before and well after the avoided crossing, thus leading 
to delocalized eigenfunctions 11 , 1^ ■ 



In conclusion, we have developed an spectroscopic 
technique to enable experimental investigation of mode- 
dynamics evolution along the chaotic transition in open 
chaotic billiards. The observed mode-dynamics evolution 
shows that openness tends to suppress avoided crossings 
compared to the closed billiard cases. We could directly 
measure the coupling strengths induced by ray chaos 
among encountering modes. Our measurements would 
serve as a valuable asset for anticipated but currently- 
nonexisting semiclassical theory for coupling strengths 
between modes in a mixed phase space. 
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